Limits on spin-dependent WIMP-nucleon interactions 
from the Cryogenic Dark Matter Search 
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The Cryogenic Dark Matter Search (CDMS) is an experiment to detect weakly interacting massive 
particles (WIMPs), which may constitute the universe's dark matter, based on their interactions with 
Ge and Si nuclei. We report the results of an analysis of data from the first two runs of CDMS at the 
Soudan Underground Laboratory in terms of spin-dependent WIMP-nucleon interactions on 73 Ge 
and 29 Si. These data exclude new regions of WIMP parameter space, including regions relevant 
to spin-dependent interpretations of the annual modulation signal reported by the DAMA/Nal 
experiment. 

PACS numbers: 95.35.+d, 14.80.Ly 



The nature of the dark matter which dominates struc- 
ture formation in our universe is one of the most press- 
ing questions of modern cosmology 0, 0, El • A promis- 
ing class of candidates is weakly interacting massive par- 
ticles (WIMPs) T], particularly the lightest neutralino 
in supersymmetric (SUSY) extensions to the Standard 
Model Q- Many groups have sought to detect WIMPs 
directly via their elastic scattering off atomic nuclei 

The nucleon coupling of a slow-moving Majorana neu- 
tralino (or of any WIMP in the extreme non-relativistic 
limit |(|) is characterized by two terms: spin-independent 
(e.g. scalar) and spin-dependent (e.g. axial vector). 
When coherence across the nucleus is taken into ac- 
count 0, these two terms behave very differently. The 
neutralino has similar scalar couplings to the proton and 
neutron 3] , and nucleon contributions interfere construc- 
tively to enhance the WIMP-nucleus elastic cross section. 
Thus, though neutralino-nucleon cross sections for such 
interactions are generally orders of magnitude smaller 
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than in the axial case |8|], scalar couplings dominate di- 
rect detection event rates in most SUSY models for ex- 
periments using heavy target nuclides. 

In contrast, the axial couplings of nucleons with op- 
posing spins interfere destructively, leaving WIMP scat- 
tering amplitudes determined roughly by the unpaired 
nucleons (if any) in the target nucleus. Spin-dependent 
WIMP couplings to nuclei thus do not benefit from a 
significant coherent enhancement, and sensitivity to such 
interactions requires the use of target nuclides with un- 
paired neutrons or protons. Spin-dependent interactions 
may nonetheless dominate direct-detection event rates in 
spin-sensitive experiments in regions of parameter space 
where the scalar coupling is strongly suppressed. This 
can provide a lower bound on the total WIMP-nucleus 
elastic cross section, since spin-dependent amplitudes are 
more robust against fine cancellations [9j- In general, 
consideration of such couplings when interpreting exper- 
imental results more fully constrains WIMP parameter 
space and allows exploration of alternative interpreta- 
tions of possible signals 0, In this work we ex- 
plore the implications of recent results from the Cryo- 
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genie Dark Matter Search (CDMS) experiment for spin- 
dependent elastic scatterin g. S ome limits with a previous 
data set have appeared in and constraints on spin- 
independent interactions with this data set are discussed 



with ao — a p + a n and a\ 



The functions 



acpci 

0. 



The Cryogenic Dark Matter Search 0, ^| seeks to 
detect WIMPs via their interaction with nuclei in semi- 
conductor crystals. CDMS uses ZIP detectors |l2( to 
discriminate between electron recoils (induced by most 
backgrounds) and nuclear recoils (induced by WIMPs 
and neutrons) on an event-by-event basis via a simulta- 
neous measurement of ionization and athermal phonons. 
Under standard assumptions about the galactic halo (de- 
scribed in Q), CDMS currently sets the strictest upper 
limits on spin-independent WIMP interactions 0, [lj| . 

The CDMS detectors are made of natural Ge or 
Si, both composed predominantly of spinless isotopes 
with negligible sensitivity to spin-dependent interactions. 
However, each contains one significant isotope with non- 
zero nuclear spin: 73 Ge (spin-9/2) makes up 7.73% of 
natural Ge, while 29 Si (spin-1/2) makes up 4.68% of nat- 
ural Si. Each isotope contains a single unpaired neutron, 
making CDMS much more sensitive to spin-dependent 
interactions with neutrons than with protons. 

The analysis presented here is based on the combina- 
tion of two data runs taken at the Soudan Underground 
Laboratory, a deep installation which provides a rock 
overburden of 2090 meters water equivalent. The first 
Soudan run used a single tower of 6 ZIP detectors (4 Ge, 
2 Si) and recoil energy thresholds of 10-20 keV. From 
October 2003 through January 2004, 52.6 live days of 
WIMP-search data were acquired. We work with the re- 
sults of the "current" analysis detailed in [l^, in which 
one candidate event was identified in Ge. The second 
run added a second tower (2 Ge, 4 Si) and acquired 74.5 
live days of data between March and August of 2004. We 
work with the 7-keV threshold analysis described in ^| i 
also with one candidate event in Ge. Both candidates are 
consistent with expected backgrounds, and no candidates 
were observed in Si. Scaling the exposures before analysis 
cuts by the isotopic abundances given above, we obtain 
a total of 11.5 (1.7) raw kg-days 73 Ge ( 29 Si) exposure. 

The spin-dependent interactions of a given WIMP are 
characterized by its WIMP-proton and WIMP-neutron 
spin-dependent couplings, a p and a n . Unlike in the spin- 
independent case, the relative strengths of these cou- 
plings may vary significantly with neutalino composition 
|24| , though t hey are often of similar magnitude in mod- 
els of interest [25| . The differential cross section for spin- 
dependent WIMP-nucleus elastic scattering at momen- 
tum transfer q can be written as |26| 
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where v is the incident WIMP velocity, J is the nuclear 
spin, Gf is Fermi's constant, and the "spin structure 
function" S{q) is given by 

S(q) = alSoo(q) + a a 1 S 01 (q) + a?Sn(g), (2) 



Sij(q) encompass the magnitude of the spin associated 
with the nucleon populations, as well as the effects of the 
spatial distribution of that spin at non-zero momentum- 
transfers. These must be determined separately for each 
nuclide using a nuclear structure model. Such models 
may be compared based on their q = nucleon spin ex- 
pectation values ({S p ) and {S n }) and the accuracy of their 
predictions for the nuclear magnetic moment. These 
models also show that spin correlations (e.g. polarization 
of the even nucleon group by the odd nucleon group) give 
a nucleus with an unpaired neutron a residual sensitivity 
to a p . 

For 29 Si, the major efforts to determine nuclear spin 
structure have been large-basis shell model simulations 
by Ressell et al. 0] and Divari et al. [2]}. The results 
of both calculations agree in the zero-momcntum-transfcr 
limit and reproduce the experimental magnetic moment 
(fi = — 0.555/ijv) to within 10%. We follow the former. 

The most complete shell model studies of the high- 
spin 73 Ge nuclide have been carried out by Ressell et al. 
[lH and Dimitrov et al. [28|. The former result requires 
"quenching" to bring its predicted value of the nuclear 
magnetic moment (/i = — 1.239/iAr) in line with experi- 
ment (/i = — 0.879^ijv), while the hybrid model used in 
the latter does not (/i = — 0.920/ijv)- Both models give 
values of (S n ) within ~ 2% of one another, but their 
values for (S p ) differ by a factor of 3. We follow Dim- 
itrov et al., but also compute Ge limits following Ressell 
et al. (the "alternate form factor") to give an indication 
of nuclear model uncertainties. Note that the structure- 
function fits given in the latter are invalid for recoil en- 
ergies En > 50 keV. We thus assume no sensitivity to 
such recoils when using this model (thereby limiting our 
sensitivity for high WIMP masses in this case). 

We follow a "model-independent" framework described 
in [Til I24I I29I ] for the interpretation of experimental re- 
sults in terms of spin-dependent interactions. For easy 
comparisons between experiments, we report allowed re- 
gions in two M x — a SD planes: one in the limit a n = 0, 
one for a p — 0. For this purpose, we express limits in 
terms of WIMP-nucleon cross sections 
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where fi xp , n is the WIMP-nucleon reduced mass. We also 
plot allowed regions in the a p — a n plane for two choices 
of WIMP mass, M x . 



Figure H shows upper- limit contours in the M x — a SD 
plane in the limiting cases of pure neutron coupling 
[cbp — 0) and pure proton coupling (a n = 0) for the 
CDMS data sets, computed using Yellin's Optimum In- 
terval method |3(j . Dashed curves in Fig. ^ and [21 show 
the limits obtained using the alternate form factor. As 
benchmarks, we also include recent limits from other 
leading experiments, as well as 3er allowed regions based 
on DAMA/Nal's reported annual modulation (lCf. The 
latter is computed following [ll| , based on the 2-4 keVee 
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FIG. 1: Upper limit contours (90% confidence level) for recent CDMS data sets, plotted in the cases of (left) pure neutron and 
(right) pure proton coupling. We show limits based on Ge (solid) and Si (dash-dot) from the combined Soudan data (black). 
Dashed curves represent Ge limits using the alternate form factor from |15|. As benchmarks, we also include interpretations 
of the DAMA/Nal annual modulation signal [lOj (filled regions are 3cr-allowed) and limits from other leading experiments: 
CRESST I [rJ as computed in ("x"s), PICASSO (squares), NAIAD [r| (circles), ZEPLIN 10 (triangles), and 
Super-Kamiokande |2(3| (asterisks; an indirect search, based on different assumptions). EDELWEISS [2l| and SIMPLE |2^| 
report limits (not shown for clarity) comparable to CDMS Si and PICASSO, respectively. As a theoretical benchmark (see 
text), horizontal dotted lines indicate expected cross sections for a heavy Majorana neutrino. Plots courtesy of [23j| . 



modulation amplitude alone and a tripling of its quoted 
la error bars. Including constraints from other energy 
bins favors somewhat lower cross sections at high WIMP 
masses, but does not modify our conclusions. 

Due to its isotopic composition, CDMS is sensitive pri- 
marily to WIMP-neutron spin-dependent couplings. The 
combined CDMS data exclude new regions of parame- 
ter space in the case of purely WIMP-neutron cou pling 
(i.e. for a p = 0). In combination with CRESST I [Iff, 
these data are inconsistent with an interpretation of the 
DAMA/Nal annual modulation amplitude in terms of 
such interactions within the standard halo model (see 
also ^2). Despite its lack of unpaired protons, CDMS 
also possesses competitive sensitivity to a p . CDMS does 
not currently set the strongest limits in the pure WIMP- 
proton case, but has begun to explore the region of this 
parameter space associated with the DAMA/Nal signal. 

To explore more general models, these results can also 
be expressed in the a p — a n plane for various choices 
of WIMP mass. Two such choices (15GeV/c 2 and 
50GeV/c 2 ) are shown in Fig. [21 The former illustrates 
the advantage of using two active isotopes: the region al- 
lowed by both nuclides (approximately the overlap of the 
corresponding ellipses) may be significantly smaller than 
either ellipse individually. In the current case, the gen- 
erally weaker limit from 29 Si serves to cut off the regions 
at large a 2 + a 2 allowed by Ge. This is significant since 
the lengths of the major axes of these ellipses depend on 
near-cancellations in the structure functions, and so have 



substantial uncertainties. Other experiments (see [Tl|) 
set more stringent constraints on a p and further reduce 
the overall allowed region. 

Specific WIMP model frameworks yield constraints on 
the relationship between a p and a n . In particular, SUSY 
neutralinos are expected to have comparable couplings 
to protons and neutrons. As an example, Bednyakov [25j 
finds limits of 0.55 < < 0.8 in the effMSSM frame- 
work, corresponding to the thin filled wedges in Fig. [21 
Imposing this constraint eliminates the overlap between 
the CDMS and DAMA/Nal allowed regions at WIMP 
masses > 25 GeV/c 2 , though compatible regions at lower 
masses remain. 

The upper limits set here do not yet constrain SUSY 
models significantly. An increase in exposure by at least 
two orders of magnitude is required to explore the most 
accessible mSUGRA models |31|. which may be possible 
with the next generation of direct detection experiments. 

One benchmark for spin-sensitive dark matter searches 
is a heavy Majorana neutrino, which has a purely axial 
coupling to nucleons. Horizontal dotted lines in Fig.^h 1- 
dicate the expected spin-dependent cross sections of such 
a particle • These lines suggest the magnitude of spin- 
dependent cross section expected from a Majorana dark 
matter candidate with weak interactions (regardless of 
production mechanism), and thereby indicate an approx- 
imate upper bound to "interesting" WIMP parameter 
space. We should note that this particle is intended as a 
benchmark rather than a serious dark matter candidate. 
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FIG. 2: Regions in the a p — a n plane allowed (at the 90% confidence level) by CDMS data. Each data set excludes the exterior 
of the corresponding ellipse. Two choices of WIMP mass M x are shown: 15GeV/c 2 {left) and 50GeV/c 2 {right). Dot-dashed 
ellipses represent Si limits, solid ellipses represent Ge limits, and dashed ellipses represent Ge limits using the alternate form 
factor. Also shown are the corresponding interpretations of the DAMA/Nal modulation signal (near-horizontal light (pink) 
filled bands are 3er-allowed). The thin dark (blue) filled wedges correspond to models satisfying 0.55 < | — | < 0.8, a constraint 
from the effMSSM framework 25]. 



Such a neutrino could be unstable, and relic density cal- 
culations suggest that it cannot be a major component 
of the dark matter if produced thermally [13 ■ CDMS 
II will begin to probe this parameter space in the near 
future, and some has already been reached by indirect 
searches. 
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